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' Variational (Rayleigh-Ritz) methods are applied to local quantum field 

(N ■ 

' of a superposition of Gaussians and the expectation value of the Hamilto- 

Q>^ \ nian is expressed in a form that can be minimized numerically. A scheme of 

successive refinements of the superposition is proposed that may converge 
to the exact functional. As an illustration, a simple numerical approxima- 
Q \ tion for the effective potential is worked out based on minimization with 

'■^ respect to five variational parameters. A variational principle is formulated 

for the fermion vacuum energy as a functional of the scalar fields to which 
the fermions are coupled. The discussion in this paper is given for scalar 
5^ \ and fermion interactions in 1+1 dimensions. The extension to higher di- 

mensions encounters a more involved structure of ultraviolet divergences 
and is deferred to future work. 
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l.INTRODUCTION 

The variational or Rayleigh-Ritz method is an important nonperturbative tool in 
quantum mechanics. In its simplest form (see e.g. ref. Q it approximates the ground 
state energy Eq of a Hamiltonian H by the expectation value 

< ^-li^l ^' > / < ^ I I- > (1) 

where ^ is a trial wave function depending on a set of parameters. The minimum of 
E\i, under variation of these parameters is the approximation to Eq. 

In general, one does not know how close to Eq the variational approximation is. It 
all depends on guessing the form of the trial function largely on the basis of intuition. 
It seems like an impossible task for systems with many or an infinite number of degrees 
of freedom like a quantum field theory (QFT). 

An obvious difficulty with QFT's is that the functional integrations required to 
compute the matrix elements in (1) are impossible unless ^ belongs to a very limited 
class of functionals essentially that of Gaussians-times-polynomials in the field variables. 
Moreover, with a polynomial multiplying the Gaussian, one obtains for each of the 
matrix elements a sum of terms proportional to different powers of the volume resulting 
in a trivial thermodynamic limit. 

Still another obstacle are the ultraviolet (UV) divergences {Eq is UV-divergent 
even for free field theories) which require the introduction of a momentum cutoff and 
of appropriate cutoff-dependent counterterms in the Hamiltonian. One then needs 
to compute the divergent parts of the counterterms exactly - for instance the exact 
divergent part of Eq must be found, otherwise the calculated value will be off by an 
infinite amount. 

Early applications of the variational approach have been largely based on a 
Gaussian wave functional. Of particular interest have been calculations of the effective 
potential ^eff{v) defined as the minimum of the expectation value of H under the 
constraint that the expectation value of the relevant scalar fiald be v. 

Various ideas have been proposed which go beyond the simple Gaussian wave func- 
tional such as the use of an appropriately unitarily transformed Gaussian Q and a 
two-particle-point-irreducible loop expansion of an effective action for local composite 
operators 

A different but, nevertheless, also promising approach is variational perturbation 
theory Q, which consists in making the separation of the Hamiltonian into free and 
interaction parts depend on a set of parameters and improving the perturbation expan- 
sion by demanding stationarity with respect to these parameters up to a given order. 
However, variational perturbation theory is not really based on a Ray leigh- type prin- 
ciple, so e.g. the effective potential calculated by it, cannot be claimed to be an upper 
bound to the exact effective potential.This is true also for a related approach employed 
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by the authors of ref. 

The present work is another attempt to go beyond the simple Gaussian functional. 
Its basic plan is to parametrize the wave functional by a set of parameters as a super- 
position of Gaussians and then to minimize the expectation value of the Hamiltonian 
numerically with respect to those parameters. 

Provided the issues mentioned above, namely the thermodynamic limit and the UV 
divergences, are successfully met, present computing capabilities warrant the hope that 
minimization with respect to a sufficiently large number of parameters can be managed 
and will produce useful results. 

In this paper we explore these possibilities in the framework of a local relativistic 
theory in 1+1 dimensions involving a real scalar field (/'(x) and a fermion field ^(x) 
with Lagrangian density 

L = L^ + L^ (2) 

= ^ - m - rip)ij (4) 

The generalization to several fields would present no new problems. However, the 
extension of the approach presented on this paper to higher dimensions encounters a 
more involved structure of ultraviolet divergences and will be taken up in future work. 

In addition to the effective potential, one other quantity for which one may hope 
to obtain a reasonable approximation in terms of a variationally determined vacuum 
functional ^ is the two-point function 

< ^'|(^(X1)(/P(X2)|* > (5) 

for all space time separations xi — X2 since from its values at equal times one can, in 
principle, obtain the weight function of its Kallen-Lehmann representation . 
This is a rather long paper. The following outline may be useful to the reader. 

1. Introduction. 

2. The form of the wave functional is discussed for a purely scalar theory. 

3. As a first orientation to the problem in its simplest form the effective potential 
for the theory is calculated by means of a single-Gaussian wave functional. 

4. A specific form for the superposition of Gaussians is chosen and the field inte- 
grations are carried out. 

5. The thermodynamic limit is taken after space has been divided into cells of equal 
size assumed to be approximately uncor related. A sequence of increasingly refined 
approximations is formulated which would arguably converge to the exact result. 
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6. Removal of the UV divergences of the scalar theory. 

7. The approach is illustrated by approximating the effective potential for the Xtp^ 
theory through numerical minimization with respect to five parameters. 

8. Variational approximations are explored for the fermion vacuum energy as a 
functional of the scalar field (p to which the fermions are coupled. 

9. Removal of UV divergences from the fermion vacuum energy. 
10. Conclusion and prospects. 



4 



3.THE FORM OF THE WAVE FUNCTIONAL. 

We consider the theory defined by the Lagrangian density in 1+1 dimensions. 
Initially, we enclose space in a box of length V with periodic boundary conditions on 
the real field ip{x). The thermodynamic limit V^oo will be taken in Section 5. 

The Schrodinger picture Hamiltonian is 

where m„ is the unrenormalized mass and A is the unrenormalized coupling constant. 
We seek to minimize the expression 

^ ^ 1 <4MM1^ (7, 

V-^oo V < > 

We take the trial functional in the form of a superposition of Gaussians 

*(^) = I3p(5)exp{-^ I dxdy{(p{x) - g{x))K{x,y){ip{y) - g{y))} (8) 

9 

where K{x, y) is a real, symmetric, positive kernel and the summation runs over some 
set of functions to be specified. Ideally, this should be a complete set in some sense. 

The virtue of form (8) is that it allows us to carry out the integrations over 
explicitly. 

Superposing Gaussians with different kernels would be incompatible with the ther- 
modynamic limit. Indeed,let Ki and K2 be two such kernels. Then < ^ \Hv \ ^ > and 
< ^'1^' > , as a result of the cp integrations would contain terms proportional to the 
determinants of {Ki + K2)^ , Ki~i and K2'^ . Assuming that the kernels Ki{x,y) and 
K2{x,y) depend only on x-y (for translation invariance), these determinants may be 
expressed, for large V, as 

V fdk, ~ ^ V fdk 



"^P^" 2" / ^log(^i + ^2)},exp{-- J ^log(ifi)} (9) 

V f dk 

exp{--y-log(K2) (10) 

where Ki and K2 are the Fourier tranforms of Ki and K2. But then the requirement 

that ratios of terms in < * \ Hv\ * > and < ^'|^' > be finite in the limit V^oo leads 

to _ _ 

fdk ,Ki+K2. _ 

which implies Ki = K2, since they are positive kernels. 

Similarly, a translationally invariant wave functional consisting of a Gaussian times 
a polynomial consisting of several terms of the form 



/ 



dxidx2dxz ■ . . F{xi — X2, X2 — X3,.. .)(p{xi)(p{x2) ■ ■ ■ (12) 
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of different degrees in (p would lead to a sum of terms proportional to different powers 
of V. 

At any rate the single-kernel superposition (8) may be adequate. One may argue 
that any "reasonable" functional, by analogy to ordinary functions,can be represented 
as a superposition of Gaussians having the same quadratic term in the exponent. 



6 



3.THE SINGLE-GAUSSIAN FUNCTIONAL 

It is instructive to consider first in some detail the case in which the wave functional 
is a single Gaussian term. We consider the wave functional 

'^{(p,v) = exp{-^ J dxdy{(p{x) - v)K{x - y){(p{y) - v)} (13) 

in which translation invariance is explicit and the average value of the field is arranged 
to be V. The expectation value of the energy density is 

= ^D-\0) + ^{-^ + ml)D{x)U=o + lmy + 

24 4 V ; g w \ J 

where 

i>(x-rf. <^lfJ;WI^> -,^ = A-V-rt (15) 
The variational equation 

(16) 



5D{x) 

can be solved explicitly. We introduce the Fourier transform 



D{k) = j dxe-'^'^D^x) (17) 



in terms of which (16) reads 

1 1 / , 9 9 X A ^ 9 r dk 



8D{k) 2 
Solving for D{k) we obtain 



{k-' + ml) + ^{v^ + l^b(k)} = (18) 



where the mass parameter (the "renormalized mass") is the solution to the "gap 
equation" 

2 9 . A o A /" dk 



< = <+2^' + 2 7 (20) 

The integral over D{k) diverges at large momenta k, so a cutoff A must be introduced. 
For ml to be finite, must depend on A so as to cancel the divergent part of the 
integral in (18). With no loss of generality we set 

The finite mass M replaces as one of the basic constants of the theory, the other 
one being A. The gap equation becomes 

1 1 /Vf2 

»;=M^ + -A.^ + -Alog(^) (22) 



"V 
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The minimum of is 

The sum of the first two terms can be shown (by calculating the contribution of tlie 
relevant perturbation diagrams) to contain the exact UV-divergent part of the energy 
density and is independent ii v. It is interesting to compare this divergent part to that 
of the one-loop approximation to Ei^f which is contained in the zero-point energy 

E^,Uoop = ll ^^k^ + ml + ... (24) 

Clearly, the variational energy is "infinitely" lower. 
We may define the finite expression 

Veff^'^\X,M,v) = E^i,{X,M,v,A) - E^{X,M,0,A) (25) 

as the approximate effective potential. It is an even function of v and it increases 
monotonically with v"^ as long as 

^ < min( ^ ^ ) = 2.160 (26) 



87rM2 X '21ogx' 

For A/SttM^ > 2.160 the approximate effective potential develops a minimum at some 
nonzero value Vq of v'^. This minimum value becomes negative for A/SttM^ > 2.439 
indicating vacuum states that break the ip —ip symmetry of the Hamiltonian. 

It is worth pointing out that for A/SvrM^ > 2.439 (symmetry breaking lower 
value for the energy density is achieved by superposing just two Gaussians, one centered 
at (f = vq and the other at </? = —vq: 

■^{(f) = cos/3*+ -Fsin/3*_ (27) 

Note that in the V ^ limit the matrix elements 

< >, < *+|(^|*- >, < > (28) 

vanish, because they are proportional to 

exp{-^^;^ J dxdyK{x-y)}. (29) 



It follows that as (3 ranges from to 7r/2 the expectation value of (p varies from vq to 
—Vq while E-in remains constant and equal to V^^^^{±.vq). Note that this "improved" 
effective potential is a convex function of as expected for the exact effective potential. 
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4.SUPERP0SITI0N 

In this Section we discuss the general superposition (8). We make a special choice 
for the set of g functions and we obtain a convenient expression for the expectation 
value so that the limit V ^ oo can be taken in the next Section. 

We begin by writing the kernel K as 

K = (A^t)-i (30) 

where ^ is a nonsingular real integral kernel to be specified. 

We make the summation over g in Eq.(8) concrete by expressing g(x) in terms of a 
real constant v and a set of parameters ci, C2, ... 

g{x) = v + ^ CaAha{x) (31) 

a 

where {ha{x)} is an orthonormal set of functions in jC2{—oo,oo). We have set 

Aha{x) = J dyA{x,y)ha{y) (32) 

Ideally, the set {ha{x)} should be complete. Realistically,though, wc shall have to 
assume that even with an incomplete set of such basis functions good results will be 
achieved provided their linear combinations can adequately represent functions that 
are localized anywhere in (configuration) space and also lie, in Fourier space, within 
the range of wave numbers characteristic of the field-theoretic system being discussed. 

The summation over g in Eq.(8) will be realized as an integration over the c pa- 
rameters 

Y.p{g)^ f{dc)p{c) (33) 

a 

where we introduced the shorthand notation 

(dc) = dReCadlmCce (34) 

a 

p{c) = p(ci,C2,...) (35) 

With respect to the set of functions {ha{x)} the essential choice is the subspace 
they span, since within the same subspace any change of basis is equivalent to a linear 
unitary transformation of the set of c variables. 

Carrying out the (p integrations wc obtain 

= {det KrU{dc){dc')p{c)p*{c')eM-lj2(^a-c*y} (36) 
<^'|^|^'> = idetK)--2 [{dc)idc')p{c)p*{c')eM~lT.ica-c^f} 

Y^Ah0{x)C-^-^) (37) 

/3 ^ 
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Note that the c variables are "deeoupled" in tlie exponential of these expressions 
and the same is true for analogous expressions for the expectation value of the product 
of any number of fields. This is the reason we have expanded g{x) in terms of the set 
{Aha} rather than {ha}- 

In what follows we shall use the shorthand notation 

^^'''^ /(dc)(dc')p(c)p*(c')exp{-|E.(ca-C,*)2} ^ > 

The condition 

< ^'^(^^'^ > = V (39) 
< *|* > 

imposes on p(c) the constraint 

<ca + c':>=0 (40) 

We introduce the following real, symmetric two-, three- and four-point correlation 
coefficients 

4^1 = <(^— )a. > (41) 

Jaia2a3 = < ( ^ )ai ( 7i )a2( 7, )«3 ^ (^2) 
- ( 



and the real symmetric kernels 

Q^^\x,y) = ^ ha,{x)4^lhaM (44) 

aia2 

We also introduce the connected four point correlation coefficient 

jc =7 _ /(+) /(+) _ /(+) /(+) _ /(+) /(+) (45) 

It vanishes whenever a subset of its indices is associated with c variables uncorrelated 
to the c variables associated with the remaining indices. Note that due to Eq.(40) J^a2 
and Jaia2a3 are already connected in this sense. 

The expectation values needed to calculate Ex^, can now be expressed in terms of 
these connected correlation coefficients. For the product of two fields we find 

±lk^^0l^=^^D(..,.,) (46, 

where D is the real kernel 

D = A{1 + Q^+^)A^ (47) 
We proceed to satisfy this operator relation by choosing for A the real kernel 

A = r>^(l + Q(+))-| (48) 
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The kernel D{xi,X2) represents an approximation to the (unrenorniahzed) two- 
point function at equal times. [Note that knowledge of the two-point function at equal 
times is sufficient to determine the weight function of its Kallen-Lehman[8] representa- 
tion and thus to determine it for all times] . We further find 

li^M|«IZ = y-.,i_Q(-.)^-.)(...) (49, 



4v Jceia2a3-^hoiiAhoi2-Ahce3 + 

«1«2«3 



Using Eqs.(46), (47), (48), (49) and (50) we arrive at the following expression for the 
expectation value of the energy density 

= \ly{dxdyD{x,y)\y=^ + mlD{x,x)} + ^rny + ^lyD-\x,x) + 
A{.4 + 3 J ^D{x,xr + 6v'J yD{x,x)} 
I J y{D-'[il + gW)^(l - Q(-))(l + Q(+))^ - l]}{x,x) + 

^{4u /"^ J2 Jaia2a3^ha^{x)Aha2{x)Aha3{x) + 
aia20;3 

j y Ja,a2a3aAhoc,ix)^ha2ix)Aha3ix)AhaAx)} (51) 



The expectation value is to be minimized with respect to D and p and the sub- 
space of functions spanned by the set {ha{x)}. Those are the "variational parameters". 
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5. THE THERMODYNAMIC LIMIT 

In this Section we introduce a factorized form for the weight function p(c) which 
allows us to take the limit F — >^ oo in our expression for 
Assuming translation invariance for D we set 

D{x,y)^D{x-y) (52) 

Then the x integrations in all but the last three terms in (51) are trivial and the y 
factors are removed. 

To remove the V dependence from the last three terms in (51) we rely on the 
connectedness of the correlation coefficients J^a2j Jaiaias, Jaia2a3a4 i-^- the fact that 
they vanish for uncorrected sets of indices. More concretely, we assume that p(c) 
factorizes 

p(c) = Ati(c)//2(c).../xjv(c) (53) 

where fj,i, fj,2---lJ,N depend on disjoint subsets si, S2, sat of the set of c variables. Then 
the connected correlation coefficients vanish unless all their indices belong to only one 
of the sets si, S2, s^r. Accordingly, it is convenient to label the c variables with two 
indices denoting by Cn^i, the i^-th variable belonging to the set s„. 

We proceed to divide space into N cells of equal length a {V = Na) and to make 
the special choice 

hn,v{x) = K{x - na) (54) 

associating the set s„ with the n-th cell. We make the cells equivalent by using the 
same function /j, for all the factors in the rhs of (53). 

The connectedness of the correlation coefficients implies that we may write 

'^ni2i,n2fc'2 ~ '^ni,n2 •^1/11)2 (^^) 
Jnivi ,n2i'2 ,n3i>3 — ^n\,n2^n2,n3Ji'i,V2,vz (^^) 

'^n\vi,n2i'2,nzV3,niV4 ~ ^ni,n2^n2,n3^n3,n4,Jvi,V2,i'3,i'A l^"^) 

Each of the last three integrals in (51) splits into N integrals which are shown to be 
identical by a shift of the integration variable by a multiple of a. If we now let TV and 
V go to infinity with fixed a, the volume V disappears from our expression the factors 
of N/V being replaced by 1/a: 

= l(-d^ + ml)D{x)\,=o + lmy + 
^ m / dxdyD~^{x - y)K^{x) 

[1 _ (1 + _ j(-))(i + j(+))i],,,,,/i,,(y) + 
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A 

— V 
6a 

A 
24^ 



XI Juuiy2,u3,ui J dxAhuj^{x)Ahu2ix)AK^{x)AK^{x) (58) 



This expression is to be minimized under variations of D, fi{c) and a. Since the 
lack of correlation between cells introduced by the factorized form of p{c) is an artificial 
constraint, we expect a = oo at the minimum. However, in actual numerical calcula- 
tions, only a finite number of cell basis functions can be used. In other words, fi{c) will 
be of the form 

^(c) /u(ci, C2, CK)n^^+i(5(i?eCi.)(5(/mci.) (59) 

Thus the value of a at the minimum will be finite and the approximate vacuum func- 
tional will not be translationally invariant. It will only be invariant under translations 
by multiples of a. 

Presumably, as the approximation gets refined with K (= number of basis functions) 
increasing to infinity, the cell size a will also grow to infinity. At any rate a value of a 
significantly larger than any characteristic length in the theory would signal that the 
effect of the assumed lack of correlation between cells (an "edge effect") is negligible. 
It would be evidence of a good approximation. 

An obvious choice for the set of basis functions {hiy{x)} is any complete orthonormal 
set of functions in C2{—a/2, a/2) which are defined to be zero outside the interval 
{—a/2, a/2), since then the orthogonality between /ii/^(x — nia) and hu2{x — n2a) for 
ni 7^ 712 would be trivially satisfied. An example is the familiar set of particle-in-a-box 
eigenfunctions 



[2 X 1 

K{x) = W - sin{i/7r(- + -)}, xe{-a/2,a/2) 
V a a 2 

= 0, x 4{-a/2,a/2) (60) 

(61) 

Linear superpositions of such functions or their derivatives, in general, are discontinuous 
at cell boundaries. An alternative, attractive choice are wavelets p. These are of the 
form 

ij^j{x) = 2ia^^V^(2-'- -n) (62) 
n,j = 0,±1,±2,... (63) 

where ip is specially constructed to ensure completeness and orthonormality 

dx'll)*^,j{x)tpn',j'{x) = Sn,n'Sj,j' (64) 
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The "mother wavelet" ^ can be chosen to possess continuous derivatives up to any 
given order. 

If ip is centered at x = with a width of order 1, then tpn,j is centered at a; = 
with a width of order a/2-', while its Fourier-space width is ~ 2^ /a. Thus, in order 
to use wavelets in the present context, we would replace the index by a set of two 
integers 

f {j, q} (65) 

which take values 

j = 0,1,2,3... 
q = 0,±1,±2,±3,... 
-2^-1 < q<2^'^ (66) 

We define 

K{3,q}(^) = 25a-5^(2-'(^ -n) -g) 

n = 0,±1,±2,... (67) 

Eq.(54) is satisfied because 

hn,{j,g} (x) = ho,{j,q} {x - ua) (68) 

Thus h^^(j^q) is centered at x ~ no + q/2^a with a width of the order of a/2^ 

The set of orthonormal functions defined by Eqs.(66) and (67) is not complete 
since wavelets with widths greater than a (i.e. those corresponding to j < ) are not 
included. This again is an "edge effect" expected to become negligible for large a. 
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6. REMOVAL OF UV DIVERGENCES 

In this Section we separate out of the additive UV-divergent term. We also 
determine the cutoff dependence of the unrenormahzed mass rriy, which renders the 
theory finite. 

The variational equation 

=0 (69) 



SD{x) 

cannot be solved explicitly for D{x) in the general case. However, its large k behaviour 
can be obtained provided the Fourier transforms of the basis functions hi,{k) vanish 
fast enough. Assuming that all terms containing hi^{k) drop out we obtain for D[k) 
the same asymptotic behaviour as that of the free theory: 

D{k)k-^±oo ^ ^ (70) 

This behaviour is indicated for the two-point function in 1+1 and 2+1 dimensions in 
perturbation theory but not in 3+1 dimensions for an interacting theory. 
Going back to we note that the quantity 

1 dE^ 1 , A 9 X f dk 



2v dv 2 " 12 4 

A 



Jui,i'2,u3 j dxAhvj^(x)Ahv^{x)Ahi,^{x) (71) 



12av 

must be finite. Therefore,assuming that the last term is UV finite, the divergent part 
of the integral over D (the "tadpole" term) must be cancelled by that of the unrenor- 
mahzed mass term. Just as in the single-Gaussian case, with no loss of generality, we 
set ^ 

thereby introducing the finite mass parameter M to take the place of m^. 

To proceed, we need to extend Eq.(70) and assume that D has an asymptotic 
expansion of the type 

D{k)k^±oo = ^(1 + ^ + 0(1^1"'')) (73) 
a>0,/3>2 

where 71, a and /? are constants. Again this is supported by perturbation theory in 
which a = 2. As a consequence of Eq.(73) we find that the expressions 

00 _ ]^ 

-00 2^^^^^^ ~ 2Vk^ + M^^ ^^^^ 

°° rlh ^ 1 ~ 1 

( JtTT + ^^'^ W - i 1^1) (75) 



/- 



-00 2Tr'8D{k) 2 ' ' 2 
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are UV finite. Tlius our expression for Eq, takes tlie form 

dk,,, 1 >^ dk 1 



= —\k\-—(M^-- y + 

' 27r' ' 2A^ 2 7_A 27r 2\/FTm2 



/ 



dk , 1 1,9 - . 1 I, IN A?;^ 
:^(^- + -^k^Dik) - - \k\) + — + 

+terms containing 

The sum of the first two terms of this expression contains the exact UV-divergent 
part of the vacuum energy density of the theory. These terms are independent of v and 
of the variational parameters D{x), n{c), a, and {h^}- We may simply drop them and 
retain the following finite expression to be minimized. 

J 2Tr^8D{k) 2 ^ ^ 2 ' 

2A^ 2 7 27r^ ^ ^ I^WTW^' 



^ X] / dxdyD ^{x-y)K^{x) 

[1 _ (1 + _ + + 

A /■ 

—V ^ Jyuv2,v3 J dxAhvi{x)Ahy^{x)Ahu:i{x) + 

H '^ui,u2,u3,ui j dxAh^^{x)Ah^^(x)Ah^^{x)Ahi,^{x) (77) 



A 



In higher-dimensional theories the assumption that the last term in the rhs of 
Eq.(71) is UV-finite cannot be maintained. This is clear even in the super-renorma- 
lizable (2+l)-dimensional case in which the UV-divergent part of m\ is calculable. In 
that model, in addition to the "tadpole" term \D{0)/2 associated with the diagram (a) 
of Fig.l, there is also the divergent contribution of the two-loop diagram (b) of Fig.l 
which must be related to the last term in Eq.(71). 
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(a) (b) 

FIG.l Diagrams contributing to the UV divergence 
of the two-point function in 2+1 dimensions. 



In 3+1 dimensions a host of interrelated UV problems must be faced. The large- 
momentum behaviour of the connected two-point function D(k) is not the same as that 
of the free theory but is determined by the anomalous dimensions of the field (provided 
an UV fixed point exists); the wave function renormalization constant is UV-divcrgcnt 
and so is the unrenormalized coupling A. It seems evident that in these theories, in 
order to achieve the minimum of Eq,, the correlation coefficients must be such as to 
make the sums in Eqs.(71) and (77) diverge in the A — oo limit. 
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7.A NUMERICAL CALCULATION OF THE EFFECTIVE POTENTIAL 

In this section we calculate numerically an approximation to the effective potential 

using the simplest superposition of Gaussians formed with just one basis function per 

cell: 

g{x) = V + ^ Cnh{x — na) (78) 

n 

and with real coefficients Cn- In this case the sum of the last three terms in Eq.(77) is 
simplified to 

^[1 - (1 - j(+))(i - J ^\h{k)\^/b{k) + 

^vJz I dx[Ah{x)f + J4 / dx[Ah{x)f (79) 

Oct J Zi^tt J 

where we have set 



= <{^?> (80) 

= <C-^f> (81) 

J4 = <{^f>-^<{^f>' (82) 
A{x-y) = L>^/2(x-y)/(l + JW)V2 (83) 

The function h{x) must satisfy the orthonormality condition 

J dxh{x — na)h{x — n'a) = 5nn> (84) 

which may be expressed in terms of the Fourier transform of h as 

I E IMfc + ^)P-l (85) 

s=—oo 

Subject to this constraint, h{x) is itself a variational parameter. However, for 
simplicity, we make a definite choice. Note that if we take the Fourier tranform of h{x) 
to be of the form _ 

^^^^ ^ (aEr=^ool/(^ + ^)P)^/' ^^^^ 
then Eq.(85) is satisfied for any f for which the sum converges (certain integrability 
conditions are obviously also necessary). For f(x) we take the second-order cardinal 
spline and find 

Mfc) = [^sin(^)]2[l-^sin2(^)]-V2 (87) 
which is the scaling function for the associated wavelet. 
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Since D{x) approximates the (unrcnormalized, equal time) two-point function it is 
appropriate to adopt the KaUen-Lehman form for it i.e. 



with the positive weight function a to be determined. For instance, approximating a 
by a sum of delta function terms, we would write 

In order to keep as few variational parameters as possible, we simply take only one 
such term 

Finally, we adopt a Gaussian-times-polynomial form for the weight function //(c) 
so that the integrations over the c's can be explicitly done: 

/x(c) = (1 + uic + U2C^ + U3C^) exp(-7cV2) (91) 

The constraint (40) can be used to express e.g. ^3 in terms of ui,U2 and 7. Thus, 
for given values of A, M and v, the expectation value of the energy density is an explicit 
function of the five variational parameters 77^,0,7,^1,^2. 

To find the minimum of a FORTRAN minimization program was run on a PC 

with a Pentium 133 CPU. Each numerical evaluation of Ei^, took .5 sec on the average. 
Each search for the minimum Eg of Eij, took ^ 2min. The product ma ranged roughly 
from 0.6 to 1.3, so the values found for Eg cannot be claimed to be close to those 
of the exact energy density. Nevertheless, they provide a rigorous upper bound to the 
corresponding values of the exact Veff{v). 

In Fig.2 the lowering Es — Eg of the vacuum energy density achieved by the su- 
perposition is plotted vs. the coupling A for t> = 0. Being zero at A = it decreases 
steadily with increasing v. 

In Fig.3 Es and Eq are plotted vs. v for A/M^ = 50 and A/M^ = 100. For 
A/M^ = 100 Es displays a minimum at v = ~ .8. 

Just as in the single-Gaussian case of Section 3, if a linear combination of "^{vq) 
and ^(—vq) were used as a wave functional the energy density in the interval {—vo,vo) 
would be lowered to the value of Eg at v = zizVQ and the curve would become convex. 
Presumably, this "extra" superposition using different values of v would not have been 
necessary, if a complete set of basis functions had been used. 
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FIG. 2 The improvement Es — Eq on the approximate vac- 
uum energy density, in units of , obtained by a su- 
perposition of Gaussian functionals vs. the couphng 
XioY V = 0. 




coupling=100 , 



0.2 0.4 0.6 0.8 



FIG. 3 Comparison of the approximate energy den- 
sities Es and Eg-, in units of M^, plotted 
against the exp. value v of the field for 
A/M2 = 50 and A/M^ = 100. 
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8.FERMI0NS 

We now consider coupling a fermion field 'ip{x), to the real scalar field ^p{x) with 
the Hamiltonian 

H = Hs + Hf (92) 

where 

/I d 
dx-lip^ix), {-ioi— +f3m + T(p{x))i;{x)] (94) 

The hermitean matrix T is some linear combination of /3 and ^/?75 with real coefficients. 
We deal explicitly with the (l+l)-dimensional theory (so we could have replaced a, 
if3'y5, and /? by the Pauli matrices (Ti,ct2 and (73) but, actually, the discussion in this 
Section is valid for any dimension (with obvious notational adjustments). 

Let Ef{ip) be the ground state energy density of Hp for a given (c-number) scalar 
field configuration (p{x). Then the ground state energy of H is the same as that of the 
Schrodinger-picture operator 

H = Hs + VEf{^) (95) 

In this Section we shall look for variational approximations to Ef{(p) by explicit func- 
tionals of (p. 

Consider a general Hamiltonian of the form 

Hf = J dxdx'^[tlj\x),h{x,x')^{x')] (96) 

where ip{x) is a fermion field operator satisfying the anticommutation relations 

{'ll)a{x),1pjj{x')} = Sa,i3S{x - x') 

{M^),M^')} = {V'i(^),V'J(^')} = o (97) 

The kernel h{x,x') may be any hermitean c-number kernel and, in particular. 



d 

h{x, x') = S{x - x'){-ia— + f3m + Tip{x)) (98) 

If {ua{x)} and {va{x)} are complete sets of positive and negative-energy eigen- 
spinors of h{x,x') and {Eq,}, {E'^} the corresponding eigenvalues, we may expand h 

as 

h = Y^iErPr + E',P;) (99) 
r 

where Pa, P^ are the projections on the corresponding eigenspinors i.e. the operators 
with kernels 

Pr{x,x') = Ur{x)ul{x') 

Prix,x') = Vr{x)vlix') (100) 
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The ground state of Hf is the state l^ij > which is annihilated by the totahty of 
the operators 

dxul{x)ip{x), / dxip^ {x)vr{x) (101) 



Assuming < >= 1 we have 

<nfmx),i;^{x')]\nf >=Y,Pr{x,x')-Y,P'Mx') (102) 

r r 

and 

<nf\Hf\nf> = ~Y,Tr{{Pr-Pl){EsPs + KPl)) 

r,s 

= ~J2(Er + mTr{l) (103) 

r 

where, in 1 + 1 dimensions, Tr{l) = 2. Prom this expression we must eventually identify 
and extract the UV-divergent parts. 

Now let |r2o > be the ground state of the free Hamiltonian 



= j dxdx'-[i)\x),h^^\x,x')tl){x')] 



h^^\x,x') = Six - x')(-ia-^ + Pm) (104) 

ox 

Clearly, Ef{(f) must be equal to the minimum value of 

^ < Wno\Hf\Wno > (105) 

for all possible choices of a unitary operator W. 

Actually, the minimum value is attained even if we limit W to the class of unitary 
operators under which the field ip{x) transforms linearly i.e. 

W^i;a{x)W = j dx'Ma,p{x,x')^p{x') (106) 

where Ma,p{x,x') is a unitary kernel. Indeed, in this case we have 

< W^lQ\Hf\WnQ ~^ S Tr{M{P^'') - P^''^')M\EsPs + KP',)) (107) 

r,s 

where , are the eigenspinor projectors of /t^^^ . 

It suffices to choose M so that it maps the subspace of positive (negative) energy 
eigenspinors of h into those of /i^^^ respectively, i.e. 

M(^pW)Mt = Y.Pr 

r r 

M{^p(^y)M^ = Y^Pr (108) 
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Then the rhs of Eq.(107) becomes equal to the ground state energy of Hf as given by 
Eq.(103). Since 

Y.P^'\x,x')-Y.P^'yM = J |expMx-xO-^|^ (109) 

we conclude that Ef{ip), the ground state energy density of Hf{if), is the minimum 
value of the expresion 

F{<p,M) = ^(|^exp(zp(x-xO^^^^^).,.Mt,(x',y) 



r,s,t,u 



d 

(-ia^ +Pm + rip{x))t,uMu,r{y, x) (110) 



under variations of the unitary kernel M. 
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9. REMOVAL OF FERMION UV DIVERGENCES 

A formal choice of M for which F{ip, M) actually attains its minimum value is 
Moller's wave matrix given by the time-ordered exponential 

.0 

S{ip) = Texp{-i dte^p{ihot)T(p{x)exp{-ihot)) (111) 

^— oo 

where x and Hq arc operators on the space of Dirac wave functions. 

Under fairly general conditions S{(p) maps eigenspinors of into eigenspinors of 
h. Therefore, according to the discussion given in the previous Section, we have 

Ef{<p) = F{^,Si<p)) (112) 




The familiar vacuum graphs of Fig.4 correspond to successive terms in the expansion 

of 

Ef{^)-Ef{0) (113) 

in powers of (p. The first term of this expansion, the term linear in is given by (with 
r = Ai/? + A2i/375) 

This is a logarithmically divergent "tadpole" contribution to Ef{ip) which can be 
cancelled by the inclusion of a c{A)(p(x) type counterterm in the Hamiltonian density. 

The next term in the expansion of Ef[(p) — Ef(0), is also logarithmically divergent. 
Its divergent piece is 

- lrr(l)(A? + Ai)i J dV(^)logA (115) 
requiring a f'^{x) counter term (mass term) in the Hamiltonian density to cancel it. 
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Unfortunately, the Moller matrix S{(p) is not available in a closed form so that, 
after we substitute F{ip,S{ip)) for Ef{(p) in Eq.(96), we can carry out the integrations 
over (p explicitly in the matrix element 

J id^)^\v)Hi<p)^i^) (116) 

with ^'((/j) a sum of Gaussian terms. We must settle for an approximation F((p, M) with 
M such that the integrations are possible. One possibility would be an exponential 
form e.g. 

M(x, x') = 6{x - x') exp{i j dyB{x, y)(p{y)} (117) 

with B{x,y) a hermitean kernel to be determined variationally. However, this would 
not reproduce the UV-divergent parts of Ef(ip) — Ef{0) exactly. Clearly, getting those 
divergent parts exactly is necessary - otherwise our approximate vacuum energy density 
would differ from the exact value by an infinite amount. 

To remedy the situation we note that these UV divergences are associated with 
the large momentum behavior of the matrix elements of T(p(x) between momentum 
eigenstates of the free Dirac hamiltonian 

<p',u'\Tip{x)\p,v>, v,v' = ±l (118) 

where \p, it > denote the positive (negative) energy eigenstates of —iot^ + jSm with 
momentum p. Thus if we define, as a modified version of Tlp{x), a cut-off kernel w{ip) 
by ^ 

<p\v'\w{'p)\p,V >= e{-\p + p'\ -Pma^) <p',v'\Tip{x)\p,V> (119) 

and construct the time-ordered exponential 

S (if) = T eyip{—i / exp(i/iot)ti)((/3) exp(— i/ioi)} (120) 

J—oo 

then F{ip,S{ip)) can be shown to contain the exact UV-divergent part of F{ip, S{ip)), 
namely 

F{ip, S{ip)) - F{ip, S{ip)) = finite > (121) 

Furthermore, by taking the momentum parameter Pmax large enough, one can make 
the expansion of F{cp, S{ip)) in powers of converge fast, so that it may be approx- 
imated by just a small number of terms (one must at least include, of course, the 
UV-divergent terms). 

In general, one cannot expect F{ip, S{p)) to be a good approximation to Ef{ip) 
(although it would still be an upper bound to it). However, one can improve on it by 
considering functionals of the form F((p, U S((p)) where U is a convenient unitary kernel 
like the simple exponential one given by Eq.(117). With U being an exponential and 
S{ip) approximated by a polynomial in f, the ip integrations in (116) could be carried 
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out explicitly. In the resulting expression the kernel B(x,y) would play the role of a 
variational parameter. 

All this can be readily extended to higher dimensions. There will be, of course, 
additional UV-divergent terms in Ef{(p) — Ef{0). For instance, in 3+1 dimensions 
there is a ip'^{x) term with a logarithmically divergent coefficient. 
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10. PROSPECTS. 

We conclude with a few remarks concerning future prospects for the variational 
approach proposed in this paper. 

At present, it appears that UV divergences are the main obstacle to extending the 
calculations described in this paper to renormalizable field theories in higher dimen- 
sions. Actually, for super-renormalizable theories in 2+1 dimensions the divergent parts 
of the counterterms are calculable since they are associated with a finite set of pertur- 
bative diagrams. In that sense they should be easier to handle. Things are considerably 
more involved for (3+l)-dimensional theories, as indicated in Section 6. 

Irrespective of whether these methods will ultimately prove applicable to realistic 
(3+l)-dimensional models (including vector fields), they could still be used to test 
other nonperturbative methods by comparing with their respective results on lower- 
dimensional systems. 

On the numerical side, a more ambitious venture than the example given in Section 
7 would require the introduction of several basis functions per cell. Let us assume, for 
example, a polynomial- times-Gaussian form for the weight function: fx{c) 

H{ci,C2, Ck) = (1 + X! ^i^^^i + X! ^ij^i^j + 

i i,j 

^flk^i^j^k + ■■■) exp{-^ ^TsCj} (122) 

The number of possible u parameters grows fast with the number of basis functions 
K. Even if we stopped at the cubic terms, there would be 55 parameters for if = 5 
and 285 parameters for K = 10. To those we must add the 7's, a and the parameters 
necessary to represent D{x). Clearly, the number of basis functions needed to achieve a 
useful result is crucial for the feasibility of the calculation; but it can only be determined 
by numerical experimentation. 

On the other hand, the virtue of any refinement in the description of the wave 
functional will be gauged by the amount by which it actually lowers the expectation 
value. This may help develop a "physical intuition" just as in variational calculations 
of atomic and molecular systems. For instance, it may turn out that the form (122) is 
not expeditious or that not all of the u coefficients in it are of equal importance. 
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